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j«^ ■ Abstract: We consider mode coupling in multimode optical fibers using 

either two Bragg gratings or a Bragg grating and a long-period grating. 

We show that the magnitude of the band edge curvature can be controlled 

Qj \ leading to a flat, quartic band-edge or to two band edges at distinct, 

nonequivalent A:-values, allowing precise control of slow light propagation. 
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Fig. 1. Schematic of tiie coupling between the forward and backward propagating modes in 
a guided-wave structure using two superimposed Bragg gratings or Bragg and long-period 
gratings. Top: schematic of mode wavenumbers. Bottom: Solid and dashed arrows indicate 
mode coupling by individual gratings. 



The ability to control the propagation of waves using the strong geometrical dispersion in 
photonic structures opens new opportunities for manipulating optical pulses. In particular, the 
regime of slow light can be achieved when the wavelength is tuned close to the edges of a 
photonic band-gap, as was demonstrated both in one-dimensional [ 1] and two-dimensional IJ] 
geometries. The properties of slow-light modes strongly depend on the shape of the dispersion 
curve: if the curvature at the band edge is eliminated, to lowest order (O — (Oe°^ k'^, where (Oe is 
the band-edge frequency and k the wavenumber This results in a flat band-edge for which the 
group velocity Vg°^ k^ °^ {(0 — (Oe)^/'^. In addition, the curvature can be made to change sign so 
that the band edges appear at distinct non-equivalent fc-values. 

A flat band edge can lead to high energy densities and to an increased density of states at 
frequencies near the band edge ElB- with associated changes in the radiative properties of 
sources. The increase in the energy density is associated with fact that for a given, low group 
velocity pulse, the range of k values in a quartic band is larger than in a quadratic band, so that 
the pulse can be more strongly compressed in space. The increased density of states directly 
follows from the flattened shape of the dispersion relation. Slow modes with non-vanishing 
phase velocities can also be used to create high-g resonators Q. It was earlier shown that such 
dispersion characteristics can be achieved in extended multi-layer anisotropic structures [S] 
|6l or omni-guide fibers HUS). Flexible control of higher-order dispersion is also possible in 
photonic-crystal waveguides |l7][8l. 

In this work, we suggest a generic approach for the engineering of flat and offset band-edges 
in conventional fibers based on the grating coupling of forward- and backward-propagating 
modes of two different symmetries, i.e., 4 modes in total. In the most common geometry in such 
structures, a Bragg grating is used to couple the fundamental mode to itself, leading to Bragg 
reflection and the opening of a (one-dimensional) photonic bandgap. The manipulation of this 
gap requires an additional degree of freedom, which is obtained by coupling two modes (such 
as the fundamental mode and a higher-order mode) using an additional grating, either a Bragg 
grating, or a long-period grating (see Fig. [T]i. This is somewhat similar to earlier proposals 
lUIll, however in that work, a thin-film geometry in which the additional field has a different 
polarization was considered. Our proposal is based on effects which cannot be achieved with a 
dual Bragg grating 121 optimized only for a single guided mode. The suggested design does not 
require the use of anisotropic layers |l3]|6l or a special fiber configuration ||4]|5|, being readily 
accessible with established grating-writing methods in conventional optical fibers. 

We consider an optical fiber supporting two different modes, .-#i and ./#2 (for example. 



LPol and LP 02)- In the vicinity of a chosen frequency Cl, their dispersion relations are given 
by ill = iV'i {Ki — Ki), Q.2 = ±V2(^2 — ^2)- Here the D.j are the frequency detunings from Q., 
Kj are corresponding wavenumber shifts, Vj are the (positive) group velocities, and Kj (> 0) 
are reference wavenumbers at Cl. The signs + and — correspond to forward and backward 
propagating waves, respectively. 

As mentioned, we introduce a grating supporting the simultaneous resonant coupling of two 
pairs of forward- and backward-propagating modes ^i and ^2- This can be achieved with a 
superimposed modulation of the refractive index along the fiber 

An{x,y,z) = AniCos[Kiz + <j)i]Rix,y)+An2Cos[K2Z + <j)2]Rix,y), (1) 

where the R{x,y) defines the photosensitive cross-section of the fiber. We choose K"i = 2Ki + 5i 
and K2= Ki+ 0X2 + ^2. Here 8j are small detunings from the two resonances, and C7 = ±1. In 
both cases a Bragg grating with wavevector K\ couples the forward- and backward-propagating 
^i mode. When a = +1 an additional Bragg grating couples the counter-propagating modes 
^\ and ^2, whereas for C7 = — 1 the co-propagating modes are coupled by a long-period 
grating. The modulation amplitudes A«i and A«2 define the coupUng strengths between the 
^\-^\ and ^\-^2 modes, respectively. Whereas it was previously found that multi-mode 
coupling and conversion with superimposed Bragg gratings can give rise to special features 
in the reflection and transmission spectra |fT0l[Tn[T2l[T3][T4l . we demonstrate below that the 
suggested superlattice modulation configuration enables new possibilities for dispersion control 
of band-edge slow-light states. 

We assume that the structure is much longer than the beat length between the two modes, 
L^ 1^2 ±^1 1 \ 1^1 1 ': 1^2! '1 which is satisfied under typical conditions. Then, the pulse 
propagation can be modeled by the coupled-mode equations ifTOl [T2l for the envelopes of the 
forward and backward propagating modes. Writing the electric field as 

E(.x,y,z,f) = {gi(.x,y)[Mi(z,f)exp(/^iz)+wi(z,f)exp(-/^iz)] 

+g2(-^,3')[M2(z,0exp(/<7A'2z) + W2(z,f)exp(-/(7A'2z)]}exp(-;i2f)+c.c., (2) 

where gi,2(-*,3') are transverse profiles for the two modes and c.c. denotes the complex conju- 
gate, we see that the forward-propagating envelopes are u\ and U2 for a = +1 or mi and W2 for 
(7 = — 1, and the backward-propagating envelopes are w\ and W2 for a = +1 or wi and 112 for 
(7 = — 1. We choose the scaling where the normalized intensities define the energy density of 
modes per unit length of the fiber, and obtain the dimensionless equations 

i^\ — V iV\ -z^ — Vp\W\ exp(/5iz + ii)\ ) + P2W2 exp(i52Z + 'fc) = 0, 
at dz 



i^r iVi -rr h PiMi exp(-;5iz - /0i) + p2M2exp(-/52Z - (^2) 

ot dz 

dii2 dui 

j^— + iC7V2^+P2Wiexp(;d2Z + i02) =0, 



(3) 



dt dz 

dw2 dw2 / c- \ ^ 

i^: ioV-)^: hp2"iexp(— ;dhz-;ffl2) = 0. 

dt dz 

Here the (real- valued) grating strengths pj are approximately given by 

Pj = -^ J J ^^yf^ix,y)giix,y) ■gj{x,y), (4) 

with A() = T-TicjCl. Note that the values of Pj can be chosen independently by se- 
lecting the strength of corresponding grating modulations and that Eqs. ^ satisfy the 



energy-conservation condition, ^(|miP + |wip + |M2p + |vv2p)/<5f = ^(^ Vi|miP + ^iKip — 
aV2\u2?- + oV2\w2\'^)/dz. 

The dispersion relation of the system can be calculated by analyzing the eigenmodes of 
Eqs. (O in the form, 

u\ = f/iexp[i(A:+5i/2)z — /(Br + /0i/2] , 

wi=Wi exp [i{k - 5i /2)z - im - iipi /2] , 

"2 = U2 exp [i{k —81/2 + 82)1 — i(Ot — i(j)i /2 + (^2] , 

W2 = W2 exp [i{k + 81/2- 82)1 - im + /^i /2 - 1(^2] ■ (5) 

We present the eigenmodes as sums of symmetric and antisymmetric vectors, Uj = Fj+ + Fj^ 
and Wj = Fj+ — Fj^ for y = 1,2. After substituting these expressions into Eqs. (O, and obtaining 
the eigenmode equations, we find that they can be formulated in the matrix form 

where 

/ («±pi)/yi-5i/2, ±p2/Vi \ 

"^'^ y ±P2CJ/V2, coa/V2 + 8i/2-82 )■ ^'^ 

Therefore, k^ are eigenvalues of a square 2x2 real-valued matrix M = M+M^. Thus there 
appear four branches of dispersion curves (Oj{k), with symmetry (aj{k) ~ (Oj{—k). For con- 
venience, we label the branches such that cOi < (O2 < (O3 < (O4. From the form of the matri- 
ces M±, we observe that the dispersion does not depend on the grating phases 0^, and iden- 
tify the symmetry property C0j{k;8i,82,Pi,P2,Vi,V2) = PiCOj{kVi/pi;0,82Vi/pi - 8iVi[l + 
aVi/V2]/[2pi],l,p2/pi,l,V2/Vi) + 8\Vi/2. Then, in the analysis and figures below we take 
5i = 0, Vi = 1 and pi = 1, since for other values the dispersion curves can be obtained by 
simple scaling. 

Away from the grating resonances, C0j{k -^ ±°o) ~ ztVi 2^, i.e., there are two pairs of 
branches with positive and negative frequency detunings. Since all (Oj{k) are real for real k, 
branches cannot merge or disappear. However, for some frequencies all k can become imag- 
inary, such that the propagation of linear waves is suppressed due to Bragg reflection. Such 
a photonic band-gap appears between the second and third branches and due to the symme- 
try of the dispersion curves, the band-gap edges appear at pairs of dispersion points labeled 
(±A:2,a)2 ) (lower) and (±A:3,(»3 ) (upper edge). The wavenumbers ±k2 and ±A:3 define the 
corresponding phase velocity of slow-light modes. 

To determine the features of band-edge dispersion, we first analyze the point with A: = 0. 
Then, Eqs. (|6]l decouple, and admit solutions in the form of purely symmetric or antisymmetric 
modes when Det{M±[cOj{k = 0)]} = 0, the solutions of which satisfy 

{COj ± 1 ) {(Oj/V2 - 082) = PI/V2. (8) 

The dispersion points (i>j{k = 0) at j = 2 or 7 = 3 correspond to the edges of the photonic band- 
gap if the second eigenvalue of M[(i)j{k — 0)] is negative, i.e. if no other propagating modes 
exist in the spectral region around (i)j{k = 0). This condition is satisfied if 

^j = {(Oj/V2 - a 82)^ + {co]-l)- 2apl/V2 < 0. (9) 

Therefore, CO- ~ COj{k = 0) and kj = if rjj < 0. By choosing the grating strength and detun- 
ing, this can be realized at either of the band-edges [see examples in Figs.lUa) and (c)], or at 
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Fig. 2. Dispersion relation for structures based on (a), (b) two Bragg gratings (cT = +1) 
and (c), (d) Bragg and long-period gratings (cj = — 1). Shading marks the band-gap, and 
filled circles indicate the band edges. For all the plots, Vi = 1, V2 = 0.95, pi = 1, P2 = 0.5, 
5i = 0. Grating detunings are (a) §2 = 0, (b) ^2 = -1.765, (c) & = 0, (d) ^2 = 0.995. 



a selected band-edge [upper band-edge in Figs.|2jb) and (d)]. In contrast, if rjj > 0, the sign 
of the curvature changes sign and the edges of the band-gap appear at points with non-zero 
wavenumbers ±A;y, as at the lower gap edges in Figs.|2b) and (d). 

A key characteristic of slow-light at the band-edges is the second-order dispersion D = 
d^co/dk^. We find that the absolute value of D[(X>j{k = 0)] is proportional to \rij\. Therefore, the 
dispersion at the gap-edge becomes quartic when rjj = as (Oj = (Oj{k = 0) +0 x k^ + 0{k'^), 
and this happens precisely at the transition of dispersion curves from the single band-edge with 
kj = to double band-edges with ±kj ^ 0. The quartic dispersion occurs at the upper gap-edges 
inFigs.|2l;b),(d). 

In Fig. [3] we present the dependence of the wavenumbers {kj) and second-order dispersion 
at the band-gap edges. We note that Eqs. (|6]l and (|7]i possess a general symmetry, (o{k\52) =■ 
—(o{k; —52), allowing us to represent simultaneously the characteristics for both the upper and 
lower band-edges. Results of calculations for superimposed Bragg gratings [Figs. |3ja),(b)] or 
Bragg and long-period gratings [Figs. [3jc),(d)] demonstrate that both the wave-number and 
dispersion of band-edge slow light can be engineered by choosing the grating detuning (52) 
and strength (P2), the values of which can be freely selected by specifying the required grating 
period and the refractive index contrast. Consistent with the theoretical predictions, the absolute 
value of D is reduced to zero [visible as white stripes in Figs.[3jb) and (d)] along the path in the 
&i— P2 plane where band-edge wavenumbers kj start to deviate from zero. 
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Fig. 3. Tuning characteristics of slow-liglrt modes for (a),(b) coupled Bragg gratings 
(C7 = 1) and (c),(d) coupled Bragg and long-period gratings (c7 = — 1). Shown are the depen- 
dencies of the absolute values of the band-edge (a,c) wavenumber and (b,d) group-velocity 
dispersion on the grating strength (P2) and detuning (+52 at upper or —82 at lower gap 
edges). Marked points A-D correspond to the dispersion plots in Figs. |2ta)-(d), respec- 
tively, where the subscripts u and / indicate the upper and lower gap edges (points A and 
C do not have subscripts as they coincide for both band-edges). Note that in (a) and (c), 
the band-edge wavenumbers are exactly zero throughout the dark blue region. Normalized 
parameters are the same as in Fig.|2] 



We note that in Eq. ([T]) we assumed that the grating planes are orthogonal to the fiber axes. 
It was pointed out by Erdogan and Sipe ifTOl that by tilting these planes the coupling between 
selected modes can be enhanced. However, the detailed design of the grating is outside the 
scope of this work. 

We stress that the technology for writing gratings in optical fibres is mature and the fabrica- 
tion of the gratings described here should not be a problem, in principle. An additional degree 
of freedom that can be brought to bear in possible experiments is that the frequency difference 
between the resonances associated with the two gratings can be tuned be stretching the fiber. 

In conclusion, we have demonstrated that both the phase velocity and the strength of 
second-order dispersion in the slow-light regime can be controlled in optical fibers with super- 
structured Bragg gratings or Bragg and long-period gratings designed for coupling of funda- 
mental and higher-order modes. In particular, it is possible to realize flat band-edge dispersion 
of quartic type, when the second-order dispersion is suppressed. As a final point, though we 
considered here an optical fiber geometry, any multi-moded guided-wave geometry in which 
gratings can be fabricated is suitable for dispersion control of slow-light. 

This work has been supported by the Australian Research Council through the Centre of 
Excellence CUDOS research projects. 



